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Dilute solutions of strongly charged polymer electrolytes undergo, upon addition of multivaltent 
salt to the solutions, a phase transition from extended conformations to collapsed or bundled ones. 
Upon further addition of salt they experience a second transition, a redissolution back into extended 
conformations. This paper presents a theoretical study of the structure and properties of the phase 
diagram of these solutions. On the basis of simple phenomenological observations a schematic phase 
diagram is constructed that allows a simple and explicit determination of the direction of the tie 
lines within the coexistence region. The actual shape of the coexistence boundary is determined by 
means of a model mean free energy functional that explicitly includes the possibility of association 
of both counterions and coions to the electrolyte. It is found that it is possible to redissolve the 
electrolytes into conformations where the bare charge of the electrolyte is overcompensated by the 
counterions but, due to the associated coions, can have either sign of total effective charge. When 
colon association is possible, the redissolution approximately coincides with the reassociation of the 
coions and counterions in the bulk of the solution. 



I. INTRODUCTION 

Dilute solutions of strongly charged polymer elec- 
trolytes (polyfilectrolytes) such as sin^ and double 
stranded DNAld~Q, polystyrene sulphonatflE, and other bi- 
ological and synthetic polymeric systemsoQ exhibit phase 
and/or structural transitions upon addition of multiva- 
lent salts to the solutions. First, the addition of mul- 
tivalent salt causes a collapse of the extended structure 
of flexible polyelectrolytes or the bundling in mono- or 
multi- molecular structures of semiflexible and rigid elec- 
trolytes that segregate into a dense phase, a phenomenon 
which in the case of DNA is also know as DNA conden- 
sation. Upon further addition of salt, the electrolytes 
redissolve acquiring again extended or unbundled states. 
There are many different examples of this type of sys- 
tem and it is clear that while the strong charge of the 
electrolytes is required for the transitions, the precise ge- 
ometry of the polymers is not a crucial factor. The fact 
that the transitions are driven only by a combination of 
electrostatic interactions with entropic and steric rfFects 
has been shown by many computer simulations.utJ 

The first transition to a collapsed state is explained 
by the onset of attractive interactions between al- 
most neutral structures formed by electrolytes whose 
charge is alippst,compensated by associated multivalent 
counterions.M~E3 The redissolution transition indicates 
that as the concentration of multivalent salt increases 
the system finds more energetically favorable conforma- 
tions without packaging its structure. The model pre- 
sented in this paper emerges from two previously pro- 
posed ideas regarding the redissolution transition. In a 
previous workEJ it was argued that a determining factor 
for the redissolution is the nature of the salty environ- 
ment into which the electrolyte redissolves. In particular, 
it was emphasized that at large multivalent salt concen- 
trations the thermodynamics of the multivalent counte- 



rions and monovalent coions is affected by the Bjerrum 
associationta where the ions of the salt within the solu- 
tion become strongly correlated with each other. Thus, 
an scenario for the redissolution is that the multivalent 
ions associated with the electrolyte prefer, without leav- 
ing the polymer, to be exposed to such environment. 
On the otlier hand, it has been argued by the group of 
ShklovskiitJ that the main reason for the redissolution is 
the decrease in energy aicailable to the system by over- 
charging the elect rolytesEJ. Since the collapsed structure 
has to be effectively neutral, the electrolytes redissolve 
so that more counter-ions can associate to them. Both 
these ideas can be combined and the result provides an 
interesting picture of different possible charged states for 
the redissolved electrolyte. 

The central point of this paper to exhibit, by means 
of a simple model, that for large amounts of multivalent 
and monovalent salt added to the system, the redissolu- 
tion transition is marked by the onset of reassociation of 
the multivalent salt ions on the bulk of the solvent. In its 
redissolved state the electrolyte acquires a large number 
of multivalent counter-ions that overcompensate its bare 
charge but, also, associated with them, a large number 
of coions that form a second associated layer. Due to the 
large amount of screening of the electrostatic interaction 
near redissolution conditions, the decrease in energy due 
to the reassociation of coions with the multivalent coun- 
terions at the electrolyte surface is of magnitude similar 
(but smaller) to the energy of reassociation in the bulk 
of the solution. Thus, this second layer of associated 
coions can only exist if Bjerrum association is already 
present in the bulk of the solution. This second layer 
enhances the decrease in energy of the electrolyte due to 
the overcharging of the first layer. The presence of the 
coions reduces the net charge of the region near the elec- 
trolyte and thus decreases the Coulombic barrier to mul- 
tivalent ions overcharge. Furthermore, in such a state. 
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there is a greater number of particles in strongly corre- 
lated states thus decreasing the energy per electrolyte 
and making its redissolution more favorable. The to- 
tal charge of the near region of the electrolyte that con- 
tains all these associated charges can have either charge 
sign, but in all cases the charge from the multivalent 
ions overcompensates the bare charge of the electrolyte. 
This multiple-species association has been observed ia. 
recent computer simulations by Tanaka and Grosber^ZI 
where the electrohqte is a small colloidal sphere and hy 
Holm and Kremeiila where the electrolyte is a rigid rodta. 
The possibility of this, second associated layer has been 
known for long timeEj'Ey but its implications for these 
phase transitions has not been fully explored. Figure 1 
presents an sketch of different possible charged states of 
the near region of an electrolyte. Consider, for exam- 
ple, a section of electrolyte with total bare charge —12, 
and multivalent counter-ions of charge z = -1-3. At low 
amounts of multivalent-salt the total charge of the sys- 
tem is negative, but partially compensated by (a) positive 
monovalent or (b) multivalent ions. The collapsed state 
(c) is neutral with the electrolyte charge compensated 
by positive multivalent ions, where different electrolytes 
or sections of electrolytes share the positive charges. In 
the redissolved state the charge of the electrolyte is over- 
compensated by multivalent ions but the total associated 
charge can be (d) positive, (e) neutral, or (f ) negative due 
to the association of a second layer of negative ions. 
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FIG. 1. Sketch of different charge states of a rod-like 
electrolyte with associated counterions and colons. At low 
amounts of multivalent-salt the total charge of the system is 
negative, but partially compensated by (a) positive monova- 
lent or (b) multivalent ions. The collapsed state(c) is neutral 
with the electrolyte charge compensated by positive multiva- 
lent ions. In the redissolved state, the total charge can be (d) 
positive, (e) neutral, or (f) negative. 

A second goal of this article is to present a simple con- 



struction of the phase diagram of the multi-species elec- 
trolyte system. While conditions for collapse or redis- 
solution implicitly assume some form of phase transition 
its global structure has not been previously considered in 
detail. This article proposes, on the basis of simple phe- 
nomcnological observations, a construction of the struc- 
ture of the phase diagram, the location of the coexisting 
phases, and the slopes of tie-lines. 

In the rest of the paper, all free energies are measured 
in units of fcT, with k the Boltzmann constant and T 
the temperature. All concentrations are presented as 
number densities. Lengths are measured with respect 
to a microscopic length unit a that is the average size of 
monomers of the polymer and ions, which are assumed 
to be of roughly the same size. The strength of the elec- 
trostatic interaction is quantified by the Bjerrum number 
B = /eakT which is energy of interaction of two ele- 
mentary charges at its nearest possible distance in water 
(or the medium considered). 



II. STRUCTURE OF THE PHASE DIAGRAM. 
A. Boundaries of the phase diagram 

The basic components of the system are: a solvent, 
a dilute polyelectrolyte with monomer concentration cj), 
monovalent salt s and a multivalent salt p. Since each of 
these components dissociate into different ions, we actu- 
ally start with six different ionic components in solution. 
To reduce this number, it is possible to consider the case 
in which the charges of the ions are —1 for each monomer 
of the polymer, +1 for the original counterions of the 
polymer, —1 and -1-1 for the ions of the monovalent salt, 
and +z and —1 for the ions from the multivalent salt. 
It is further assumed that all of the +1 ions are identi- 
cal, and similarly for the —1 ions. Thus there are only 
four components, besides the solvent: dissociated poly- 
mer at concentration p — (j), positive multivalent ions 
m — p, positive monovalent ions t — -I- s, and nega- 
tive monovalent ions u — zm + s. A further constraint is 
that each phase formed has to be electroneutral so that 
the condition —p + zm + t — u = is always satisfied. 
These restrictions allow for a description of all the ionic 
species in terms of only three concentrations, either of 
salts (4>,p,s), or ions {p,ra,t). 

The allowed values of the concentrations of ions can 
be described as follows. Each ionic particle density 
should be non-negative, and in the reduced representa- 
tion (p, m, t) the requirement m > becomes the con- 
straint t — p -\- zm > 0. Carrying this constraint to the 
salt representation shows that ((/>, p, s) should satisfy 



O<0, 
0<p, 
min(zp, (j)) < s. 



(1) 
(2) 
(3) 
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This last inequality has the following interpretation. 
Consider the case of a state with a small amount of sol- 
vent and a dense mixture of only polymer ions with mul- 
tivalent ions, so that (p, m, t) = (x, — x/z, 0) (and u — 0). 
This state is physically realizable but its description in 
terms of the salt variables is (0, p, s) — {x, x/z, ~x). To 
obtain this state one can start with a mixture of polyelec- 
trolyte and multivalent salt and then extract monovalent 
salt, that is, extract both positive monovalent counteri- 
ons (given by the polyelectrolyte) , and negative mono- 
valent ions (given by the salt). The boundaries of the 
possible states of the system in the salt variables repre- 
sentation are sketched in Figure 2. 




t 

FIG. 2. These diagrams show sections of the planes that 
bound the region of allowed states for the system in the 
(p,m,t), and (</>, p, s) representations. As explained in the 
text, the allowed states extend into negative values of the 
monovalent salt variable s. 



B. Location of the dense phase 

Experiments have shown the phase segregation of neu- 
tral complexes with little amount of solvent from the bulk 
of the dilute polyelectrolyte solution. The precise loca- 
tion of this dense phase in the salt or ionic components di- 
agram has not been precisely determined experimentally. 
Here it is only necessary for our purposes to assume that 
this phase has zero or very small amount of monovalent 
ions present and that is rich in polyelectrolyte neutralized 
with multivalent ions with a very small amount of sol- 
vent. The dense phase is not a unique point in the phase 
diagram, but it can be described as a region in the neigh- 
borhood of the point (p, TO, t) — ($, <i>/z, 0),or {(j),p,s) — 
(<i>, —<i>), where $, is the (large) concentration of 
monomers in this phase and, as explained above, the neg- 
ative value of s reflects the depletion of monovalent ions 
from these states. To simplify the model, however, it 
can be assumed that the dense phase consists of only one 
point, namely {4>, p, s) = (<&, — $). 



C. The tie lines 

The precipitation of the collapsed neutral structures 
from the dilute solution indicates that the system has 
entered into a coexisting phase. A point (</>, p, s) in- 
side the coexistence region decomposes into two coexist- 
ing points with compositions (<&, <l>/z, — $), namely the 



dense state, and a second polymer dilute point of the 
form {(j) — x^, p — x^/z, s -\- x^). The assumption that 
most of the high polymer density points of the region of 
coexistence lie in a very small region and are modeled as 
one single point immediately gives rise to this specific de- 
composition, and it is clear that the slope of the tie- lines 
is always determined by this construction. 

To summarize, the phenomenological observation that 
the dense polymer phase is neutral with the polyelec- 
trolyte charge neutralized by multivalent salts leads to 
a phase diagram of the form shown in Figure 3 includ- 
ing the tic lines shown. (The choice of cross section for 
different values of the t variable are meant to aid visual- 
ization of the diagram and do not reflect actual energetic 
considerations. An example of actual cross-sections in 
logarithmic variables is presented in the Figure 4). 
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FIG. 3. Scheme of the region of coexistence in the reduced 
description of the components of the system. It is assumed 
that all tie lines connect points of the bounding surface of the 
coexistence region with a single point ($,"I>/2:,0) at the one 
extreme of the coexistence region. Four tie lines are shown for 
the points corresponding to a section defined by a constant 
value of the t component. 



III. FREE ENERGY FUNCTIONAL 

To be able to write down an effective free energy func- 
tional that describes the polyelectrolyte phase transition 
it is necessary to consider coherently but separately the 
contributions to the free energy from three different re- 
gions. There is a dense precipitated phase that is rela- 
tively poor in solvent and rich in polymer and multivalent 
ions. Within the dilute phase, there are two distinct re- 
gions: the one occupied by the polymer and its associated 
charges, and the region away from the polymer occupied 
only by solvent and ions. 

There are many variants of the problem all of which 
present similar behavior. The same type of phase dia- 
gram occurs for flexible and rigid polyelectrolytes of dif- 
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ferent lengths. It is simplest to write a model for a col- 
lection of rigid rods (say short or medium size DNA). In 
this case the main barrier arising from the structure of 
the polyelectrolyte to the formation of a collapsed state 
is purely entropic. If one considers very long semiflexible 
molecules, the collapse is monomolecular and the barrier 
is of a elastic nature instead of entropic, and requires 
the introduction of further parameters. The results are 
nevertheless similar. 



A. Two States modeling 

While it is clear that the region of phase coexistence 
does not separate two topologically disconnected regions 
in the full phase diagram, it is nevertheless useful to write 
down two separate free energy functionals for the states 
at the ends of a tie line. It is very hard to write down one 
single free energy functional that interpolates between 
the dense and dilute phase, especially for flexible elec- 
trolytes. When using different functionals for each phase, 
a phase transition can be properly determined only if 
both functionals are good approximations of the absolute 
free energy and refer to a same (arbitrary) zero energy 
point. In this case this goal is achieved by using similar 
approximations in the calculation of the free energy in 
the two states. 



B. The Dense Phase 

In the dense phase, the charges are very close to each 
other although it is likely that they retain a large amount 
of water. This phase is electroneutral and effectively 
forms a bulk region where the all the charges are strongly: 
correlated with their near neighbors. In previous workEj 
it was argued that a suitable free energy per unit volume 
Fd for this state, is given by: 

Frf($,$/z,-$) = -M$. (4) 

where M in an effective Madelung constant that depends 
on the specific geometry of the polymer and ions used but 
that reflects the strongly correlated state of the charges 
in that phase. 

C. The dilute region 

The dilute phase is not homogeneous at a microscopic 
level. The charge distribution is very different near the 
polyelectrolyte compared to regions away from it. The 
free energy of this phase can be approximated by sep- 
arating a region near the polymer from the rest of the 
solvent. Since the polymer concentration is very di- 
lute, this region can be considered as a separate ther- 
modynamic system in equilibrium with an ion reservoir. 
The amount of charge per monomer associated to the 



polymer from each of the ionic species is fm, fu, ft- 
The average ionic concentrations in the bulk of the so- 
lution excluding regions occupied by polymer is then 

(to, t, u)e// = (to - fmP/z, t - ftP, U - fuP)- 

The free energy density for the system polyelectrolyte 
plus solvent can be written as: 

Ft^^\n^+c^F,+Fs, (5) 

where Ft is the total free energy density in the dilute 
phase, Fe is the energy per electrolyte monomer due to 
the interactions between its own and associated charges, 
and Fs is the energy of the solvent an the ions in solution. 
The values of the associated fractions fi are determined 
by minimization of this expression with respect to them. 

The salty environment can be modeled as an ionic so- 
lution obeying a Debye-Huckel free energy where it is 
necessary to consider not only the multivalent positive, 
monovalent positive and negative ions but also an extra 
species, the reassociated ions. At the typical concentra- 
tion at which the redissolution of the electrolyte takes 
place, there are important corrections to the free energy 
of the solution due to this phenomenon. The simplest 
way to describe the association is to assume that a fi- 
nite fraction of all the multivalent ions in solution are 
associated with the same number of positive monovalent 
ions thus forming a new species of reduced charge (z — 1) 
with concentration TOq = /tTOg//. While in general there 
can be further reassociation to states of charge z — 2 etc., 
only the first reassociation will be treated explicitly here. 
Upon reassociation, the magnitude of the energy of a pair 
of associated ions is gt- The new concentrations of non 
associated ions are TO„ = mef f —ma, and u„ = u^f f—nia- 
The free energy for the solvent is then: 

-Fs = V'i In V'i - -^K^ - gb^a (6) 

^ — ' iZTT 

where k is the screening length of the system and is given 
by 

A«2 = ^47rz2^„ (7) 

i 

where the ionic densities ipi with charges Zi are to^, to„, 
teff and Un- 

The free energy of the near region of the polyelec- 
trolyte, per monomer of electrolyte can be expressed as 
an expansion up to quadratic order on the fractions of 
associated charges. A good choice for this functional is 

G2 

Fe = —q — gmfm — 91 ft — Qxifu + 9cfmfu- (8) 

To argue for such functional it is convenient to separate 
within the near region of the electrolyte, the contribu- 
tions to the free energy coming from a small region near 
each positive multivalent or monovalent ion, and the in- 
teraction of between these ions with others further away. 
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The long range screened interaction between charges 
gives rise to the quadratic contribution G'2g^/2, where 
the effective charge per monomer ia q — 1 — /m + fu — ft- 
The coefficient of interaction depends on the geometry of 
the system and for simphcity is chosen to have the form 
G2 = GiB ln((l + k)/k), where Gi is a geometric factor. 
For large screening lengths ^ 1, this term reflects 
the interaction of each charge with the logarithmic field 
of a charged rod of length k~^. At very small screening 
lengths < 1 the function goes to 0. 

When the effective electrolyte charge is near zero the 
linear terms in the free energy reflect the correlation en- 
ergy of each charge with its near neighbors. At zero effec- 
tive charge, the amount of energy gained from the addi- 
tion of a unit charge from a given ion is simply —gmi —gs, 
or —gt- The association of a positive monovalent charge 
has as prerequisite the presence of multivalent negative 
ions associated to the electrolyte, so that /„ < /,„. The 
reduction in energy for this association is required to be 
comparable (but smaller) to the energy of association for 
these pairs of ions in the bulk so that gu < gb the positive 
charge is brought into contact with a positive associated 
ion but it is also in the neighborhood of the negatively 
charged polyelectrolyte. 

The reassociation of the multivalent ions and negative 
colons influences the strong correlations between associ- 
ated multivalent ions. As a larger number of monovalent 
ions enter the near region, the multivalent ions reduce 
its correlations since their repulsion is diminished by the 
presence of the negative ions associated with them. This 
element is reflected in the correlation reduction term gc- 
More such terms describing corrections to the correla- 
tions can be written but it turns out that the inclusion 
of just this one is sufficient to recover a phase diagram 
that is similar to those observed in experiments. Omis- 
sion of this term leads to much more extended regions of 
coexistence. 



IV. DETERMINATION OF COEXISTENCE 
REGIONS 

A. Equilibrium between phases 

A state with concentrations {p, m, t) decomposes into 
two states if the following free energy, for some value of 
cc, < cc < 1, is smaller than the total energy of that 
state. 

Fcomp = xFd + Ft{<i) - a;$, p - x^/z, s + < i^t(</), p, s). 

(9) 

The boundary of the coexistence region is determined by 
points at which the condition above is an equality and 
the required value of x is zero. Expansion of Fcomp with 
respect to x leads to 

-M -^i<0 (10) 



where the effective chemical potential pre monomer for 
the addition of a neutral electrolyte complex into the di- 
lute region is 

ti = d4,Ft + il/z)dpFt~d,Ft. (11) 

On the basis of the assumption made, the final condi- 
tion for the determination of the coexistence boundary 
does not explicitly depend on the concentrations in the 
dense phase given by <i>, but only on the bulk energy per 
electrolyte charge — M and the slope of the tie-line. The 
specific value of <& is, nevertheless, important to deter- 
mine the condensed electrolyte fraction x. 



B. Collapse 



In a typical experiment performed to study the phase 
transitions of this type of system, the amount of multi- 
valent salt present is changed from to some maximum 
value in small increments. The first transition occurs (for 
low monovalent salt) when the amount of charge provided 
by the multivalent salt is approximately the same as the 
amount of polyelectrolyte. If there is a large amount of 
monovalent salt, the transition requires a larger amount 
of multivalent salt. Thus, it is possible to approximate 
the location of the transition as follows: the condition 
Eq.^ is satisfied almost as soon as full charge compen- 
sation of the electrolyte by multivalent ions occur. The 
chemical potential /i™ of the ions in the solution, disre- 
garding the Debye Huckel contribution, is ln(m — /mi^) 
with /m ~ 1- Thus, for low monovalent salt concentra- 
tions the condition is simply that there must be enough 
ions to compensate the bare charge of the electrolyte: 

zm^cj). (12) 

At larger monovalent salt concentrations, the multivalent 
ions must first displace associated positive monovalent 
ions from the electrolyte. The reduction in the energy 
from the association of the monovalent and multivalent 
ions is, per charge, approximately the same and thus the 
chemical potential in the solvent for the multivalent ions 
should be smaller than that of the monovalent ions. The 
condition for the charging of the electrolyte with mul- 
tivalent ions, and a first approximation to the collapse 
transition condition is: 

-lnTO«ln(s + (/)). (13) 
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C. Redissolution 



In solutions where even at large salt concentrations the 
free energy of the ions is dominated by their translational 
entropy and there is little correlation between charges, 
the redissolution transition occurs when the decrease in 
energy per electrolyte due to multivalent ions overcharge 
is larger than the bulk correlation energy of the dense 
phase. However, if many monovalent ions are present, 
the onset of association of the positive monovalent with 
negative multivalent ions marks a sharp decrease in the 
energy of the electrolyte. Not only there is a decrease in 
the energy of the positive charges but also it is possible 
to associate many more ions. The repulsive interaction 
is decreased as the near region can be near neutral while 
having acquired large number of positive ions. 

To first approximation, the onset of the association of 
positive ions occurs when the magnitude of the transla- 
tional entropy of the free monovalent ions is comparable 
to the energy of association gi,. At larger monovalent ion 
concentrations, the charging of the electrolyte by positive 
multivalent ions together with monovalent colons occurs 
when the energy of association gi, is capable of compen- 
sating the translational entropy of the monovalent ions. 
Thus we have 

Inu = ln(zm + s) « —gt- (14) 

This expression provides a rough approximation for the 
location of the second transition. For fixed amounts of 
added monovalent salt, the transition location is approx- 
imately a horizontal line in the {p,m) or (0, p) plane. 



D. A numerical example. 



The results of a numerical minimization of the model 
free energy appear in Figure 3. There it is shown, for 
different values of the monovalent salt concentration, the 
predicted regions of coexistence in the {(f>, p) plane. The 
parameters used for this example are B = 3.0, gm = 9.0, 
= 9.0, gt = 9.6, 56 = 1.5, gc = 1.8, M = 9.05. The 
curve containing points a, and b has the lowest amount of 
added monovalent salt s — 1.810^^ (intended to repre- 
sent a 0.01 M molar concentration), the other two curves 
are for concentrations s — 9.010"^, and s — 1.810"'* 
(0.05 M, and 0.1 M). The largest salt concentration leads 
to smaller coexistence regions. The marked pairs of 
points (a, 6), (c, d), and (e, /) correspond to the loca- 
tion of the pairs of coUapse-redissolution transitions that 
occur when for fixed amounts of monovalent salts and 
fixed amounts of polymer, multivalent salt is added to 
the system. 
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FIG. 4. The three solid curves correspond to the coexis- 
tence region boundary for three sections of the phase dia- 
gram in the {4>,p,s) presentation, for different amounts of 
monovalent salt, details in text. The largest amount of salt 
produces the smallest coexistence region. All states outside 
the boundary for each section have a negative net charge. For 
the case with largest amount of monovalent salt the dotted 
curve marks the border of the region where the multivalent 
ions overcompensate the bare charge of the electrolyte. 

In this example, at the boundary of the rcdissolved re- 
gion for the different values of monovalent salt added, the 
states d, e and / are all overcharged with positive multi- 
valent ions but the net charge of all these states is nega- 
tive. For the choice of parameters made, all states of the 
system that are not decomposed (lying outside the coex- 
istence curve for fixed amount of monovalent salt) have a 
total negative charge. Different choices of parameters can 
lead to different charged states and, for example, in sys- 
tems where the energy of reassociation of the multivalent 
salt is very small, all the dissociated states should have a 
positive effective charge, and the phase diagraipJias the 
shape proposed in work of Nguyen an ShklovskiiEil, where 
the lines of redissolution and collapse do not join to form 
a bounded coexistence region. The model presented here 
also predicts that there is a maximum value of monova- 
lent salt for which there is a coexistence region. As in the 
scheme of Figure 2, the coexistence region has a limited 
extent in the s, t ov u axes. 

For the case of largest amount of monovalent salt 
added, (with coexistence curve containing the points 
(e,/)), the dotted line divides the region of states that, 
while effectively negative, have enough multivalent ions 
to overcompensate the electrolyte bare charge. Separat- 
ing curves for other salt concentration show similar be- 
havior, always originating near the " tip" of the diagram. 
The coexistence region sections are closed regions but 
their "left" boundaries lie at very small concentrations 
that are typically experimentally inaccessible, and are 
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not shown in Figure 4. 

In extremely dilute solutions there is no large difference 
between the original concentration of ions in the bulk of 
the solution and the effective concentrations obtained by 
subtracting the associated ions. As the amount of poly- 
mer increases, the depletion from the bulk reduces the 
amount of available ions and states with similar amounts 
of associated ions appear only at higher multivalent salt 
concentrations. This leads to an upward curving trend 
in the redissolution curves. 



specific ionic species with the parameters used. Finally, 
the addition of further ingredients that can directly mod- 
ify the chemical potential of the different charged species 
should also induce dramatic changes in the location a of 
the transitions. From the point of view of the approach 
presented here, these modifications would be reflected in 
changes of the phenomenological parameters used in the 
model. 
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V. CONCLUSIONS AND FUTURE WORK 



It is still an important but unfinished problem to re- 
fine the description of the layer of associated ions by both 
providing a concrete definition of the layer and by testing 
possible (mean field) models of the free energy contribu- 
tion of the layer to the total energy of the system. In the 
context of integral equation approaches this goal can be 
achieved by locaiing specific properties of the ion distri- 
bution functionslia, but ideally such definition ought to 
be model independent and experimentally useful. Never- 
theless, even with arbitrary definitions of the layer, a free 
energy expansion of the form Eq. H is a useful instrument 
for the analysis of this type of system. 

One would like to reduce the large number of param- 
eters used in the present description of the system as 
many properties of real systems cannot be precisely mea- 
sured or have a poor experimental definition. To use the 
ideas presented here in the context of the real system re- 
quires the collapse of information about the interaction 
between molecules and monomers into a few simple con- 
stants. These quantities refiect the average radii of ions 
and are supposed to include the steric effects of water as 
well as the specific location of charges in the electrolytes. 
A careful parametrization of simulation results, on the 
basis of the charge of the near region of an electrolyte 
wEis. for example, presented in the work of Messina et 
al.cj. In model systems studied though simulations or in- 
tegral equations is possible to reduce the number of these 
parameters and calculate them on the basis of simpler ge- 
ometric parameters. For comparison with experimental 
data it is, however, preferable to retain these constants 
as adjustable parameters. 

A couple experimental tasks are suggested by this work 
as important to improve our understanding of these sys- 
tems: the presentation of the data available from the 
point of view of phase separation of ionic species (in- 
stead of salts), a more systematic mapping of the ef- 
fective charge of the electrolytes in the redissolved re- 
gion of the phase diagram, and determination of the 
ionic composition of the associated layer (in this respect 
there has beetu acpgress by means of osmotic pressure 
measurementfcJtj). Future analysis of experimental 
data using the model presented in this article should pro- 
vide a map, even if only a rough one, of the properties of 
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